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ABSTRACT 


Investigation  of  oceanographic  data  acquired  in  the  Columbia  River 
estuary,  which  lies  between  the  states  of  Washington  and  Oregon, 
revealed  significant  variation  in  sound  velocity  with  respect  to  time  and 
location.  National  Ocean  Service  (NOS)  standards  for  echo-sounder 
acquired  soundings  require  the  knowledge  of  sound  velocity  to  within 
plus  or  minus  4  meters  per  second.  To  meet  this  requirement  in  the 
Columbia  River  estuary,  methods  were  devised  to  allow  prediction  of 
sound  velocity  based  on  location  and  the  height  of  tide. 

Confidence  intervals  associated  with  sound  velocity  predictions 
provided  by  regression  analyses  showed  a  substantial  improvement  in 
accuracy  when  compared  to  using  a  single  average  velocity  for  each 
location.  Nevertheless,  there  were  still  some  cases  where  the  NOS 
requirement  could  not  be  met. 
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I.  INTRODUCTION 


A.  NOS  HYDROGRAPHY 

"Hydrography  is  that  branch  of  physical  oceanography  dealing  with 
the  measurement  and  definition  of  the  configuration  of  the  bottoms  and 
adjacent  land  areas  of  oceans,  lakes,  rivers,  harbors,  and  other  water 
forms  on  earth"  [Ref.  1:  p.  1-3].  The  National  Ocean  Service  (NOS) 
is  the  federal  agency  responsible  for  producing  accurate  nautical  charts 
of  nearshore  areas  of  the  United  States  and  its  territorial  possesions. 
The  principal  objective  of  most  hydrographic  surveys  conducted  by  NOS 
is  to  obtain  basic  data  for  the  compilation  of  nautical  charts  with 
emphasis  on  the  features  that  may  affect  safe  navigation  [Ref.  1:  p. 
1-3].  For  this  thesis  data  analysis  is  limited  to  the  study  of  sound 
velocity  variations  and  the  correction  of  echo  soundings  for  NOS 
hydrographic  surveys.  The  type  of  data  used  for  this  investigation  is 
often  the  basis  for  other  physical  studies  such  as  the  dynamics  of 
oceanic  mixing  and  salt  transport.  No  attempt  will  be  made  to  analyze 
or  explain  such  processes. 

1 .  Sound  Velocity  and  Sound  Speed 

Manuals  and  papers  about  hydrographic  surveying  often  use 
the  term  "sound  velocity"  to  designate  what  a  physicist  would  define  as 
"sound  speed."  The  term  "sound  velocity"  will  be  used  in  this  thesis 
as  many  of  the  references  also  use  the  term. 

2.  Echo  Sounding 

Echo  sounders  are  the  basic  tool  of  the  NOS  hydrographer  for 
determination  of  depth.  Echo  sounders  measure  the  time  required  for  a 
sound  wave  to  travel  from  its  point  of  origin  to  the  bottom  and  return. 
They  also  convert  this  measured  time  to  distance  or  depth  [Ref.  1:  p. 
AG-1].  The  basic  formula  for  echo  sounder  depth  determination  is; 
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Z  =  VT/2 


(1.1) 


where  Z  is  the  depth,  V  is  the  calibrated  velocity  of  sound  for  the  echo 
sounder,  and  T  is  the  time  between  transmission  and  reception  of  the 
sound  signal.  For  depths  to  be  accurate  precise  measurements  of  time 
and  a  knowledge  of  the  mean,  or  harmonic,  sound  velocity  (Chapter  II) 
are  mandatory. 

B.  THE  VELOCITY  OF  SOUND  IN  SEA  WATER 

Wilson  [Ref.  2:  p.  641]  developed  empirical  equations  based  on 
laboratory  measurements  of  .sound  velocity  under  tightly  controlled 
conditions.  The  basic  equation  is: 

V  =  1449.14  ♦  Vj  ♦  Vp  ♦  Vs  ♦  Vsjp  (1.2) 

where  V  is  the  velocity  of  sound  in  meters  per  second,  V-p  is  a 
correction  for  temperature,  Vp  is  a  correction  for  pressure.  Vs  is  a 
correction  for  salinity,  and  Vsjp  is  a  simultaneous  correction  for  the 
combined  effects  of  temperature,  salinity,  and  depth  [Ref.  1:  p.  4-75]. 
Wilson's  equations  (Chapter  II),  which  are  currently  used  by  NOS, 
were  chosen  to  analyze  the  data  for  this  thesis. 

C.  NOS  ACCURACY  REQUIREMENTS 

The  sound  velocity  at  which  NOS  echo  sounders  are  calibrated  is 
800  fathoms  per  second.  To  obtain  accurate  depths  for  nautical  charts 
the  hydrographer  must  correct  all  soundings  where  the  velocity  of 
sound  deviates  from  the  calibrated  value.  There  are  several  methods 
available  for  determining  corrections  for  soundings,  and  they  will  be 
discussed  in  the  following  section.  The  NOS  accuracy  requirements 
[Ref.  1:  p.  4-70]  are: 

For  use  in  correcting  echo  soundings,  the  velocity  of  sound  must 
be  known  with  sufficient  accuracy  to  ensure  that  no  sounding  will 
be  in  error  by  as  much  as  0.25  percent  of  the  depth  from  this 
cause  alone.  Therefore,  the  mean  velocity  must  be  known  to 


within  plus  or  minus  4  meters  per  second.  To  satisfy  these 
accuracy  requirements,  one  must  know  the  mean  temperature  of  the 
water  to  an  accuracy  of  plus  or  minus  1  degree  Celsius  and  the 
salinity  to  within  plus  or  minus  1  part  per  thousand.  A  sufficient 
number  of  salinity  and  temperature  observations  must  be  made  so 
that  the  velocity  can  be  determined  within  the  specified  accuracy 
over  the  entire  area  sounded. 


D.  PROCEDURES  FOR  CORRECTING  SOUNDINGS 

Echo  soundings  can  be  corrected  by  either  a  direct  comparison  or 
an  oceanographic  procedure.  Direct  comparisons  involve  comparing  echo 
soundings  to  known  depths;  this  does  not  require  a  knowledge  of  the 
sound  velocity.  Procedures  which  involve  the  calculation  of  sound 
velocity  from  measurements  of  water  characteristics  are  classified  as 
"^oceanographic  and  limnologic  determinations"  [Ref.  1:  p.  4-74].  For 
this  thesis  the  direct  measurement  of  sound  velocity  by  a  sound 
velocimeter  is  also  classified  as  an  oceanographic  procedure. 

1 .  Direct  Comparison  Systems  * 

Two  procedures  used  to  determine  depth  corrections  by  the 
direct  comparison  method  are  the  bar  check  and  the  vertical  cast.  In 
these  the  effects  of  sound  velocity  variation  and  instrument  error  are 
included  in  the  results  [Ref.  1:  p.  4-71]. 

a.  Bar  Check 

The  bar  check  procedure  involves  lowering  a  bar  over 
the  side  of  a  survey  vessel  and  positioning  the  bar  under  the 
transducer  at  a  known  depth  indicated  by  markings  on  the  lines 
supporting  the  bar.  Comparisons  between  the  depth  as  indicated  by 
the  echo  sounder  and  the  known  depth  are  made  at  5-  or  6-foot 
increments  as  the  bar  is  lowered  and  raised  though  the  water  column. 
Due  to  the  handling  of  the  bar,  this  procedure  is  usually  limited  to 
smaller  vessels  and  the  quality  of  the  results  is  a  function  of  sea 
conditions.  Observations  during  rough  seas  are  subject  to  unacceptable 
magnitudes  of  error  [Ref.  1:  p.  4-71]. 


b.  Vertical  Cast 


The  vertical  cast  requires  simultaneous  sounding  of  the 
bottom  using  an  echo  sounder  and  a  lead  line.  This  procedure,  which 
can  be  accomplished  by  practically  all  survey  vessels,  is  not  as 
accurate  as  the  bar  check,  but  vertical  casts  can  be  used  to  determine 
or  verify  instrument  errors.  Where  sounding  sheets  junction,  there  are 
often  discrepancies  and  soundings  fail  to  agree.  Vertical  casts  can 
sometimes  help  resolve  such  problems  [Ref.  1:  p.  4-72]. 

2.  Oceanographic  Sound  Velocity  Determinations 

The  more  traditional  oceanographic  procedures  involve  sound 
velocity  determination  by  measurement  of  the  variables  which  enter 
Equation  1.2,  namely  temperature  and  salinity  as  functions  of  depth. 

The  Nansen  cast  and  electronic  sensors  such  as  the  CTD  (conductivity, 
temperature,  and  depth)  or  the  STD  (salinity,  temperature,  and  depth) 
are  examples  of  equipment  used  for  these  procedures.  An  alternative  is 
the  direct  measurement  of  sdund  velocity  using  a  sound  velocimeter. 

The  electronic  sensors  record  data  almost  continuously  with  depth. 

[Ref.  1:  p.  4-74]. 

To  obtain  velocity  corrections  with  Nansen,  STD,  or  CTD 
data,  Umbach  [Ref.  1:  p.  4-75]  recommends  using  the  "summation  of 
layers"  method,  which  extends  from  the  concept  of  mean  or  harmonic 
sound  velocity. 

Most  of  the  electronic  sensors  which  are  currently  in  use  yield 
precise  measurements  in  excess  of  NOS  accuracy  requirements  for  sound 
velocity  corrections  if  properly  calibrated  and  maintained  [Ref.  1;  p. 
4-74].  According  to  NOS  standards  all  electronic  sensors  must  be 
periodically  checked  for  accuracy  against  Nansen  cast  or  equivalent 
observations . 

E.  SOUND  VELOCITY  VARIATION 

The  effect  of  changes  in  temperature,  salinity,  and  pressure  on 
sound  velocity  depends  upon  the  initial  conditions  of  these  variables. 
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Wilson's  equations  are  used  to  illustrate  in  Table  I  the  effect  of 
changing  variables  by  one  unit. 


TABLE  I 

Variation  of  Sound  Velocity 


PRES 

TEMP 

SAL 

V 

delta  V 

db 

g/kg 

m/s 

m/s 

1.00 

15.00 

30.00 

1502.90 

1.00 

14.00 

30.00 

1499.30 

3.60 

1.00 

5.00 

30.00 

1464.27 

1.00 

4.00 

30.00 

1460.03 

4.25 

1.00 

10.00 

35.00 

1490.76 

1.00 

10.00 

34.00 

1489 . 48 

1.28 

1.00 

10.00 

1.00 

1449.01 

1.00 

10.00 

0.00 

1447.84 

1.17 

60.00 

10.00 

30.00 

1485.33 

59.00 

10.00 

30.00 

1485.32 

0.02 

1.00 

10.00 

30.00 

1484.38 

0.00 

10.00 

30.00 

1484.36 

0.02 

Velocity  (V) 

was  computed 

using  Wilson' 

s  equatii 

Corrections  for  sound  velocity  variations  in  estuaries  can  present 
the  hydrographer  with  a  difficult  task.  According  to  Cameron  and 
Pritchard  [Ref.  3:  p.  306],  an  estuary  is  a  semi-enclosed  coastal  body 
of  water  having  a  free  connection  with  the  open  sea  and  within  which 
the  sea  water  is  measurably  diluted  with  fresh  water  deriving  from  land 
drainage.  Estuaries  are  extremely  dynamic  systems.  Topography, 
river  flow,  and  tidal  action  are  important  factors  that  influence  the  rate 
and  extent  of  mixing  salt  and  fresh  water  [Ref.  4:  p.  4].  For  the 
hydrographer,  estuaries  are  areas  where  sound  velocity  varies  not  only 


with  location  but  also  in  time  due  to  the  tidal  cycle  and  variations  in 
river  flow.  Although  no  specific  procedures  are  given,  Umbach 
[Ref.  1:  p.  4-75]  notes  the  possible  complexity  of  such  areas: 


Under  special  circumstances  such  as  ...  in  areas  of  extensive 
estuarine  discharge  of  fresh  water  .  .  .  series  of  regional  velocity 
curves  are  required.  Regional  curves  must  be  carefully  studied  to 
determine  how  they  can  be  best  grouped  by  area  or  time  or  to 
permit  drawing  of  average  curves. 


F.  OBJECTIVES 

One  objective  of  this  thesis  was  to  investigate  the  magnitude  of 
sound  velocity  variation  associated  with  significant  changes  in  water 
characteristics  within  an  estuary. 

A  second  objective  was  to  investigate  the  feasibility  of  determining 
velocity  corrections  according  to  NOS  standards  in  areas  where  extreme 
changes  in  water  characteristics  occur  with  respect  to  space  and  time. 

A  third  objective  was  to  write  a  FORTRAN  program  to  enable 
further  testing  by  the  field  hydrographer  equipped  with  a  personal 
computer.  The  use  of  commercial  subroutines  has  been  avoided. 
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II.  METHODOLOGY  AND  PROCEDURE 

A.  AREA  OF  STUDY 

The  Columbia  River  estuary  which  discharges  into  the  Pacific 
Ocean  at  latitude  46°15'N,  longitude  124®00'W  has  variations  in  water 
characteristics  which  make  it  an  excellent  area  for  investigation. 
According  to  Dyer  [Ref.  4:  p.  96],  the  Columbia  River  estuary  is 
characterized  by  an  unusual  combination  of  large  river  discharge  and 
strong  tidal  currents.  Recent  investigations  have  revealed  differences 
in  salinity  of  up  to  30  grams  per  kilogram  at  individual  stations  during 
a  single  tidal  cycle  [Ref.  5:  p.  104].  Due  to  this  large  difference  and 
the  availability  of  data  from  a  recently  completed  survey,  this  region 
was  chosen  for  investigation. 


B.  DATA  SOURCES 

To  analyze  the  problem  completely,  water  characteristics 
(temperature  and  salinity  as  functions  of  depth)  must  be  known  at 
stations  spaced  throughout  the  region,  and  each  station  occupied  during 
a  complete  lunar  tidal  cycle.  Variations  caused  by  seasonal  changes  in 
fresh-water  runoff  should  be  examined  thoroughly. 

1 .  Columbia  River  Estuary  Data  Development  Program 

The  Columbia  River  Estuary  Data  Development  Program 
(CREDDP)  was  a  5-year,  federally  funded  study  which  included 
approximately  2  years  of  data  acquisition  for  investigation  of  physical 
and  biological  characteristics.  CREDDP  was  designed  to  meet  the  needs 
of  two  major  groups.  One  group  consisted  of  local,  state,  and  federal 
government  agencies  involved  in  planning  and  permitting  activities. 

The  second  group  included  research  scientists  and  educational 
institutions  [Ref.  6:  p.  v] , 


One  accomplishment  of  CREDDP  was  the  acquisition  of  current 
velocity,  conductivity,  temperature,  and  depth  (VCTD)  data  for  a 
better  understanding  of  physical  processes  such  as  currents,  sediment 
transport,  and  salinity  intrusion.  VCTD  profiling  was  carried  out  for 
13  consecutive  days  from  17  stations  for  periods  of  from  12  to  48  hours 
per  station.  Between  October  16  and  27,  1980,  approximately  1600 
VCTD  casts  were  made.  Of  these,  458  casts  were  made  available  to 
users  [Ref.  6:  p.  L-8] . 

Data  were  supplied  to  NPS  by  the  U.S.  Army  Corps  of 
Engineers  (COE)  at  Portland,  Oregon,  on  two  9-track  magnetic  tapes 
[Ref.  6:  pp.  1-5]. 

For  an  analysis  of  sound  velocity  variation,  VCTD  casts  from 
n  of  the  17  stations  were  investigated.  The  stations,  along  with  their 
positions  and  times  are  listed  in  Table  II.  The  plot  of  the  geographical 
positions  of  the  stations  is  shown  in  Figure  2.1.  These  stations 
covered  the  area  where  the  maximum  sound  velocity  variation  would  be 
expected  as  deduced  from  salinity  intrusion  diagrams  of  the  estuary 
[Ref.  5:  p.  104]. 

2.  NOS  Bar  Check  Data 

NOS  bar  check  data  were  used  for  a  sound  velocity  variation 
analysis.  The  descriptive  report  for  survey  H-8422  listed  the  day  each 
bar  check  was  taken  and  the  correction  for  each  5-foot  increment  of 
depth.  Also  given  was  a  general  term  for  the  stage  of  the  tide  such  as 
ebb,  flood,  high  water,  or  low  water.  Data  for  46  bar  checks  were 
acquired  between  May  9  and  July  16,  1958  [Ref.  7:  p.  1]. 

3.  NOS  Predicted  Tide  Data 

NOS  predicted  tide  data  were  obtained  for  correlation  and 
regression  with  the  CREDDP  and  the  H-8422  bar  check  data.  Predicted 
tides  were  used  because  field  hydrographers  rarely  have  the  time  or 
the  proper  tools  for  a  detailed  analysis  of  real  tide  data. 
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TABLE  II 

Columbia  River  Estuary  Data  Development  Program 


STATION  (STN)  INFORMATION 


STN 

POSITION 

ACQUISITION 

START 

TIMES  (PDT^ 

2N 

46n5' 45"N 
124*03 '31''W 

10/18/80 

0930 

10/19/80 

0230 

2S 

46*15' 13"N 
124*03 '20"W 

10/18/80 

1230 

10/19/80 

0930 

4NA 

46*14' 03"N 
123*56'09"W 

10/20/80 

2115 

10/2 1/gO 

1745 

4NB 

46*14' 47"N 
123*55' 53"W 

10/20/80 

2105 

10/21/80 

1730 

4SA 

46*12'42''N 

123°57'09"W 

10/26/80 

0007 

' 10/27/80 

0832 

4SB 

46*12' 51"N 
123*56' 55"W 

10/25/80 

2350 

10/27/80 

0820 

SNA 

46*13 ' 55"N 
123*53' 17"W 

10/20/80 

1733 

10/22/80 

0930 

5NB 

46*14' 18"N 
123*53' 14"W 

10/20/80 

2035 

10/22/80 

0913 

5NC 

46°14'07"N 
123*53' 14"W 

10/21/80 

0935 

10/22/80 

1003 

5SA 

46°11'55"N 
123*48' 38"W 

10/24/80 

0830 

10/25/80 

2331 

5SB 

46°12'05"N 
123*48' 54"W 

10/24/80 

1126 

10/25/80 

1130 

10/25/80  1130 
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Figure  2.1  CREDDP  Stations. 

C.  CREDDP  DATA  ANALYSIS 
1 .  COE  VCTD  Data  Tapes 

COE  VCTD  data  tapes  contained  pertinent  information  such  as 
header  data,  time,  temperature,  depth,  conductivity,  and  salinity. 
Mercier  [Ref.  6:  p.  L-8]  shows  a  breakdown  of  the  data  components  for 
one  record  of  data.  The  tapes  contained  more  data  records  than  were 
necessary  for  this  investigation.  The  VCTD  sampling  device  recorded 
data  every  0.04  seconds,  thereby  making  the  average  number  of 


•  ■-  V.v  .v 


records  for  a  single  cast  approximately  1500.  Having  270  casts  to 
analyze,  a  method  was  devised  to  reduce  the  data  to  a  more  efficient 
size.  A  computer  program  was  formulated  to  input  the  data  for  each 
cast  and  then  output  the  first  record  and  the  records  which 
corresponded  to  the  nearest  whole  meter  of  depth.  For  example,  a 
10-meter  cast  would  theoretically  have  11  records  of  data;  the  first 
record  would  represent  the  surface  layer,  and  the  remaining  10  records 
would  be  for  1 -meter  increments.  In  accordance  with  the  NOS  policy  of 
using  bar  check  data  measured  at  5-  or  6-foot  increments,  data  were 
eliminated  except  at  1 -meter  increments.  A  sample  of  a  profile  after 
editing  is  shown  in  Table  III.  Standard  depths  which  represent  depths 
of  1 ,  3,  5,  7,  9,  11,  and  13  meters,  were  adopted  to  further  simplify 
data  processing  and  description. 

Several  drawbacks  were  encountered  using  the  CRfeODP  data. 
Because  of  the  shallow  depths  of  CREDDP  stations,  many  of  the  casts 
were  inadequate  for  correcting  soundings  in  the  deeper  areas  of  the 
estuary.  Some  casts  did  not  include  surface  layer  data  (0  to  2 
meters) . 

2.  Sound  Velocity  Computations 
a.  Wilson's  Equations 

Sound  velocity  was  computed  for  each  record  by  Wilson's 
equations  which  use  temperature,  pressure,  and  salinity  for  the 
computation.  In  the  following  equations  V  is  sound  velocity  in  meters 
per  second;  T  is  temperature  in  degrees  Celsius,  where  -4  <  T  <30 
°C;  P  is  pressure  in  kilograms  per  square  centimeter,  where  1  <  P  < 
1000  kg/cm  '  and  S  is  salinity  in  grams  per  kilogram,  where  0  <  S  <  37 
g/kg. 

V  =  1449.14  ♦  Vj  ♦  Vp  ♦Vs  ♦  VsTp  (2.1) 

Vj  =  4. 5721 T  -  4.4532x10’2t2  - 

2.6045x10'3t3  ♦  7.9851x10'®T'^ 


(2.2) 


=  1.60272x10  'P  ♦  1.0268x10  ♦ 

3.5216x10’9p3  -  3.3603xl0'^2p4 


Vg  =  1.39799(5-35)  ♦  1 .69202xl0’^(S-35)2 


VsTp  =  (S-35)(-1.1244x10'2t  ♦  7.771 1  xlO'^T^  ♦ 
7.7016xl0'5p  -  1.2943xl0'^p2  ♦ 

3.1580x10'®PT  ♦  1.5790x10'^Pt2)  ♦ 
P(-1.8607x10'^T  ♦  7.4812x10'®t2  ♦ 
4.5283x10'^t2)  ♦  p2(-2.5294x10’7t  ♦ 
1.8563x10'3t2)  ♦  p2(-1 .9646x10’'*0t) 

The  standard  deviation  for  the  velocity  is  0.30  meters  per 
second  [Ref.  8:  p.  1357], 


TABLE  III 

4SA--Cast  6--Sound  Velocity  Variables 


COLUMBIA  RIVER 

STATION  =  4SA  .  CAST  =  031B08 

LATITUDE  46°12'42"N  LONGITUDE  123°57'09"W 

DATE  OF  CAST  1980  OCT  26  TIME  OF  CAST  5:16:  0 


TIME 

PRES 

CTD 

DEPTH 

TgMP 

COND 

SAL 

s 

db 

m 

ratio 

g/kg 

0.04 

-0.04 

0.10 

12.26 

0.4531 

15.71 

5.48 

0.85 

1.00 

12.04 

0.5252 

18.57 

12.04 

1.86 

2.01 

11.79 

0.7073 

25.91 

18.48 

2.85 

3.00 

11.78 

0.7571 

27.95 

24.60 

3.86 

4.01 

11.78 

0.7932 

29.42 

31.52 

4.86 

5.00 

11.79 

0.8184 

30.46 

37.68 

5.90 

6.04 

11.79 

0.8164 

30.37 

43.80 

6.86 

7.00 

11.80 

0.8232 

30.64 

50.48 

7.88 

8.01 

11.81 

0.8291 

30.88 

56.80 

8.89 

9.01 

11.81 

0.8302 

30.92 

63.32 

9.88 

10.00 

11.81 

0.8330 

31.04 

69.92 

10.91 

11.03 

11.82 

0.8338 

31.07 

76.12 

11.90 

12.01 

11.82 

0.8356 

31.14 

82.84 

12.90 

13.01 

11.82 

0.8346 

31.10 

b.  Mean  Sound  Velocity 


Because  echo  sounders  are  time  measuring  devices, 
distance  or  depth  is  one-half  the  round  trip  travel  time  of  a  sound 
pulse  multiplied  by  a  suitable  mean  velocity.  The  mean  sound  velocity 
(MSV)  is  given  by 

MSV  =  Z/T  (2.6) 

where  Z  is  the  distance  traveled  in  time  T.  To  correct  for  sound 
velocity  variations,  the  water  column  is  divided  into  a  series  of  layers, 
each  with  thickness  AZj  and  an  associated  velocity  Vj.  Therefore,  the 
time  interval  for  the  sound  wave  to  pass  through  the  i^^  layer  is: 

ATj  =  AZj/Vj  (2.7) 

Summing  all  layers  results  in: 
n 

T=.2ATj  =.2AZj/V:  (2.8) 

1=1  '  1=1  '  ' 


or: 


T  =  /odZ/V  (2.9) 

Substituting  into  Equation  2.6  results  in: 

MSV  =  Z/(/^dZ/V)  (2.10) 

Equation  2.10  was  evaluated  using  the  trapezoidal  method  [Ref.  9:  p. 
86].  NOS  uses  mean  sound  velocities  in  the  "summation  of  layers” 
method  described  by  Umbach  [Ref.  1:  p.  4-74]  and  Wallace  [Ref.  10; 
p.  RK530.1]. 

Two  sound  velocity  listings  computed  using  the  data  in 
Table  III  are  shown  in  Table  IV.  The  sound  velocity  column  was 


computed  using  Equation  2.1,  and  the  mean  velocity  column  was 
computed  using  Equation  2.10. 


TABLE  IV 

4SA-CAST 

6--Sound  Velocity  Data 

DATE: 

26  OCTOBER  1980 

TIME: 

0520 

(PDT) 

SOUND 

MEAN 

DEPTH 

VELOCITY 

VELOCITY 

m 

m/s 

m/s 

0.10 

1475.5 

1475.5 

1.00 

1478.1 

1476.8 

2.01 

1486.1 

1479.6  [ 

•  3.00 

1488.6 

1482.2  1 

4.01 

1490.5 

1484.1  1 

5.00 

1491.8 

1485.5 

6.04 

1491.7 

1486.6 

7.00 

1492.1 

1487.3 

8.01 

1492.5 

1488.0 

9.01 

1492.5 

1488.4 

10.00 

1492 . 7 

1488.9 

11.03 

1492.8 

1489.3 

12.01 

1492.9 

1489.5 

13.01 

1492.9 

1489.8 

14.01 

1492.9 

1490.0 

_ 

To  simplify  the  analysis  of  the  variation  of  mean  sound 
velocity  as  a  function  of  time,  a  computer  program  was  written  to 
arrange  the  data  in  individual  arrays  for  each  depth  at  each  station. 

3.  Mean  Sound  Velocity  and  Tide  Height 

Inspection  of  the  mean  sound  velocity  arrays  suggested  that 
velocities  were  maximum  at  high  tide  and  minimum  at  low  tide.  To  test 
this,  a  correlation  and  regression  analysis  between  the  CREDDP  derived 
mean  sound  velocities  and  predicted  tide  heights  was  performed. 
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a.  Predicted  Tide  Heights 

To  perform  the  correlation  and  regression  analysis,  a 
method  had  to  be  developed  to  calculate  the  height  of  the  tide  at  any 
desired  time.  One  method  would  be  to  obtain  the  harmonic  constants 
for  prediction  of  the  Columbia  River  estuary  tides  from  the  NOS  Tides 
and  Water  Levels  Branch.  These  constants  would  then  be  used  in  the 
appropriate  equations  [Ref.  11:  p.  123].  From  a  computer  programming 
standpoint  this  method  would  be  the  simplest.  A  second  method  would 
be  to  formulate  a  computer  program  based  on  the  graphical  method 
(Figure  2.2)  described  in  the  NOS  Hydrographic  Manual  [Ref.  1:  p. 
4-68].  This  method  was  chosen  because  normally  the  field 
hydrographer  is  not  supplied  with  the  harmonic  constants  for  project 
areas.  However,  the  tide  tables  are  an  essential  tool  of  the 
hydrographer. 

The  present  NOS  HYDROPLOT  computer  program  for  tide 
prediction  lacks  documentation  and  is  not  designed  to  solve  this 
particular  problem  [Ref.  10:  p.  500.1].  Therefore,  a  program  was 
written. 

Points  A  and  E  (Figure  2.2)  represent  the  predicted  low 
and  high  tide,  respectively.  The  straight  line  between  these  points 
was  divided  into  four  equal  parts  by  points  B,  C,  and  D.  Point  B’  was 
plotted  directly  below  B  at  one-tenth  the  total  range,  and  point  D'  was 
plotted  directly  above  D  using  the  same  distance  [Ref.  1:  p.  4-68], 
forming  three  regions  for  the  solution  of  predicted  heights  between 
points  A  and  E.  For  the  region  between  points  A  and  B’,  the  curve 
was  assumed  to  be  parabolic,  and  in  order  to  calculate  the  tide  height 
for  a  designated  time  in  this  region,  the  equation 

Hj  =  aTj2  ♦  bTj  ♦  c  (2.11) 

was  used,  where  H  is  the  tide  height  and  T  is  the  designated  time.  To 
calculate  the  coefficients  a,  b,  and  c;  the  equations 

Ha  =  ♦  bT^  *  c  (2.12) 
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Predicted  Tide  Graph. 
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c 


(2.13) 


Hg'  ~  * 

0  =  2aT^  >  b  (2.14) 

were  used.  For  times  that  lie  between  points  D'  and  E,  the  same 
method,  was  used.  Between  points  B’  and  D'  the  change  in  height  was 
considered  to  be  linear  with  respect  to  time  and  can  be  computed  by 
the  equation 

Hj  =  Hg.  ♦  m(Tj-Tg.)  (2.15) 

where  m  is  the  slope  defined  by 

m  =  (H0'-Hg.)/(TD-Tg.)  (2.16) 

Using  the  program  based  on  the  above  equations,  tide  heights  were 
calculated  for  the  times  of  all  VCTD  casts  used  in  the  analysis.  "Data 
Array"  means  a  set  of  sound  velocity  and  predicted  tide  height  data 
that  shows  the  variation  of  both  variables  with  time  for  a  particular 
station  and  depth  (Table  V). 

b.  Data  Plots 

Computer  generated  plots  of  predicted  tide  height  and 
observed  mean  sound  velocity  (Figure  2.3)  were  made  for  all  stations 
using  mean  velocities  based  on  standard  depths. 

c.  Correlation  Analysis 

A  correlation  analysis  routine  for  mean  sound  velocity 
and  tide  height  was  run  for  all  data  arrays.  In  addition  to  computing 
the  correlation  coefficient,  the  routine  also  computed  the  averages  and 
the  standard  deviations  of  the  mean  sound  velocity  and  the  predicted 
tide  height  for  each  data  array.  Averages  were  computed  by  the 
equation 

n 

2,  x/n  (2.17) 

1  =  1  ' 
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TABLE  V 

4SB--0ata  Array- -Depth  =  7  m 


TIDE 

MEAN 

TIME 

HEIGHT 

VELOCITY 

hr 

ft 

m/s 

0.12 

1.2 

1469.6 

1.50 

4.1 

1476.6 

2.63 

6.5 

1486.5 

3.60 

7.6 

1490.8 

4.62 

7.6 

1490.8 

5.27 

7.1 

1487.3 

6.32 

5.6 

1483.1 

9.23 

1.7 

1472.1 

11.55 

3.0 

1471.0 

12.58 

4.9 

1472.2 

13.17 

6.0 

1476.9 

17.53 

7.7 

1489.6 

21.52 

-0.8 

1471.4 

22.43 

-1.5 

1468.6 

23.35 

-1.3 

1467.7 

24.28 

-0.2 

1467.1 

25.33 

1.9 

1470.0 

26.42 

4.1 

1473.6 

27.57 

6.2 

1480.0 

28.48 

7.2 

1485.9 

29.48 

7.2 

1491.9 

30.50 

6.4 

1488.6 

31.53 

5.1 

1485.4 

32.53 

3.8 

1477.3 

where  Xj  is  the  observation  and  n  is  the  number  of  observations 


[Ref.  12:  p.  124].  Standard  deviations  [Ref.  12:  p.  128]  were 
computed  using  the  formula 


/n 


(Xj-x)2/(n-l ) 


(2.18) 


Computation  of  the  simple  correlation  coefficient  [Ref.  12:  p.  258]  of 
the  mean  sound  velocity  and  predicted  tide  height  arrays  was  done  with 
equation 


r  = 


/ 


[n^.gi  XjYj  -  ( .jSi  ^j)] 


(2.19) 


-  (.|:^Xi)2][n.|:^Yi2-  (.2  Yj^)] 
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A  sample  of  the  results  of  a  correlation  analysis  is  listed  in  Table  VI . 


TABLE  VI 

4SB--Stati$tical  Oata->Depth  =  7  m 

AVERAGE  TIDE  HEIGHT  .  4.4  ft 

STANDARD  DEVIATION  OF  TIDE  HEIGHT  .  3 . 0  f t 

AVERAGE  MEAN  VELOCITY  .  1480.8  m/s 

STANDARD  DEVIATION  OF  MEAN  VELOCITY  _  8.0  m/s 

CORRELATION  COEFFICIENT  .  0.93 


d.  Regression  Analysis 

To  determine  the  feasibility  of  predicting  mean  sound 
velocity,  a  first-order  (or  linear)  regression  analysis  was  performed  on 
all  data  arrays.  The  data  were  fitted  to  the  equation 

y  =  a  ♦  bx  (2.20) 

where  y  is  the  estimated  mean  velocity,  a  is  the  coefficient,  b  is  the 
the  intercept,  and  x  is  the  tide  height.  The  coefficient  and  intercept 
estimates  were  determined  to  minimize  the  sum  of  the  residuals 


where  Yj  is  the  observed  mean  sound  velocity  and  yj  is  the  estimated 
velocity  [Ref.  13:  p.  533].  Determination  of  the  coefficients  and  the 
intercepts  was  accomplished  by  the  normal  equations 


(2.22) 


■I  ,"1  ."-"T  •  ■.■»  T»:v»v'* 


n  n 

an  ♦  b  2  Xj  =  .2,  Yj 
1=1  '  i=l 

a  2  Xj  >  bi  Xj^  =  2  XjYj  (2.23) 

where  n  is  the  number  of  observations  for  a  particular  station  and 
depth  [Ref.  13:  p.  533].  The  sum  of  the  squares  of  the  residuals  is 

SS  =  2  (Yi-yi)2  (2.24) 

The  standard  error  of  estimate  is 

see  ^(SS  /  v)  (2.25) 

where  v  is  the  degrees  of  freedom  given  by 

V  =  p  -  n  (2.26) 

4 

where  n  is  the  number  of  sample  observations  and  p  is  the  number  of 
estimated  variables  (Ref.  14;  p.  131].  For  a  first-order  analysis,  there 
are  two  estimated  variables.  Results  of  a  first-order  analysis  are  listed 
in  T able  VI I . 

Confidence  intervals  associated  with  predicted  mean 
velocities  obtained  by  equations  derived  from  the  linear  regression 
analysis  were  calculated  to  determine  the  accuracy  of  predictions. 
Confidence  interval  computation  requires  the  standard  error  of  estimate, 
the  value  of  the  student-t  (t)  distribution  which  is  based  on  the 
desired  percent  of  confidence,  and  the  degrees  of  freedom  for  the 
sample.  For  this  thesis,  90-percent  confidence  intervals  were  computed 
using  the  equation 

Cl  =  t  t(v,0.95)[1  ♦d/n)  ♦ 
n 

[(xj-x)^  /  2  (x:-x)^]]see 
i  =  l 


(2.27) 


TABLE  VII 

4SB--Lineap  Regression- -Depth  =  7  m 


COEFFICIENT  =  2.8 

Y  INTERCEPT  =  1468.0 

MEAN 
VELOCITY 
CAST  RESIDUAL 


1 

-0.3 

2 

0.2 

3 

0.2 

4 

3.2 

5 

2.8 

6 

0.1 

7 

-5.1 

8 

-0.8 

9 

-1.7 

10 

-1.2 

11 

-1.0 

12 

-1.4 

13 

0.1 

14 

-0.9 

15 

-4.5 

16 

-4.0 

17, 

1.4 

18 

5.0 

19 

7.2 

20 

5.5 

21 

2.9 

22 

0.0 

23 

-3.2 

24 

-1.2 

25 

0.0 

26 

5.1 

27 

4.8 

28 

0.0 

29 

-1.7 

30 

-0.8 

SS 

244.8 

see 

3.0 

A  RESIDUAL  IS  THE  DIFFERENCE  BETWEEN  THE  OBSERVED 
MEAN  VELOCITY  AND  THE  PREDICTED  MEAN  VELOCITY. 


where  Cl  is  the  confidence  interval  for  a  single  estimated  velocity  (yj), 
Xj  is  the  tide  height,  "x  is  the  mean  tide  height,  and  n  is  the  number  of 
sample  observations  [Ref.  15:  p.  30]. 


Several  plots  of  the  data  arrays  indicated  a  curvilinear 
relationship  between  mean  sound  velocity  and  predicted  tide  height. 
Therefore,  second-  and  third-order  polynomial  regression  analyses  were 
performed  on  all  data  arrays  to  provide  a  better  fit  of  the  data, 
thereby  increasing  the  accuracy  of  the  mean  sound  velocity  prediction. 
The  normal  equation  method  used  in  the  first-order  analysis  was 
extended  for  calculation  of  a  second-order  polynomial.  To  check  the 
results  of  the  programs  for  the  normal  equation  method,  a  subroutine 
using  the  orthogonal  polynomial  method  was  obtained  from  the  NPS 
Computer  System  Library  [Ref.  16:  p.  8].  A  third-order  regresssion 
analysis  was  also  performed  by  the  orthogonal  polynomial  subroutine  to 
investigate  the  possibility  of  increasing  the  accuracy  of  prediction  and 
fit.  Based  on  results  of  the  investigation  and  check  (Chapter  III),  the 
first-order  and  second-order  normal  equation  methods  were  designated 
as  the  primary  methods  of  analysis.  The  normal  equation  method 
appeared  to  be  the  simplest  to  program.  The  second-order  analysis 
uses  the  equation 

4 

y  =  aQ  ♦  aix  *  a2x2  (2.28) 


where  y  represents  the  estimated  mean  sound  velocity  and  x  is  the 
predicted  tide  height  [Ref.  13:  p.  535].  Coefficients  (aQ,  a^^  ,  and  a2 ) 
were  determined  to  minimize  the  sum  of  residuals  by  the  following 
normal  equations 


n  n  .,  n 

inH  ♦  a,  2  X  .  ♦  a  2  X.  =  2  Y. 

U  =  i  1  2i  =  i  1  i  =  i  1 


3  - 


X.  Y  . 
1  1 


•7  ”3  A  '  O 

an2xi*ai2x,-*a72)’<i'*  =  2x.  '^Y. 


(2.29) 

(2.30) 

(2.31) 


where  n  is  the  number  of  casts  for  a  particular  station  and  depth 
[Ref.  13:  p.  536].  Results  of  a  second-order  analysis  are  listed  in 
Table  VIII. 
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TABLE  VIII 

4SB-- Polynomial  Regress  ion --Depth  =  7  m 


COEFFICIENTS 


an 

=  1470.0 

aV 

=  1.6 

a2 

=  0.1 

MEAN 

VELOCITY 

CAST 

RESIDUAL 

1 

-1.7 

2 

0.8 

3 

1.0 

4 

3.2 

5 

2.7 

6 

0.2 

7 

-4X7 

8 

0.6 

9 

-1.0 

10 

-1.7 

11 

-1.4 

12 

-0.9 

13 

-0.7 

14 

-2.4 

15 

-5.0 

16 

-3.2 

17. 

2.0 

18 

4.2 

19 

5.3 

20 

1.4 

21 

-0.6 

22 

-1.7 

23 

-3.0 

24 

-0.4 

25 

0.8 

26 

5.4 

27 

5.0 

28 

0.7 

29 

-0.9 

30 

0.0 

SS 

209.8 

see 

2.8 

A  RESIDUAL  IS  THE  DIFFERENCE  BETWEEN  THE  OBSERVED 
MEAN  VELOCITY  AND  THE  PREDICTED  MEAN  VELOCITY. 


Confidence  intervals  for  predicted  velocities  were 
computed  using  Equation  2.27.  The  degrees  of  freedom  for  a 
second-order  analysis  was  less  than  a  first-order  analysis  because  of 
the  addition  of  one  estimated  variable. 


D.  NOS  BAR  CHECK  DATA 


The  NOS  bar  check  data  consisted  of  46  bar  checks  which'  were 
acquired  to  determine  the  corrections  for  soundings  of  survey  H-8422. 
The  format  of  the  data  was  a  listing  of  corrections  for  5-foot  increments 
of  depth.  The  initial  analysis  included  determination  of  the  maximum, 
minimum,  and  mean  correction  for  each  depth  (Table  IX). 


TABLE  IX 

Bar  Check  Correction  Data  For  Survey  H-8422 


OBS 

MEAN 

ft 

UPPER 

LIMIT 

ft 

LOWER 

LIMIT 

ft 

RANGE 

ft 

DEPTH 

ft 

88 

-1.0 

-0.6 

-1.4 

0.8 

5 

90 

-0.8 

-0.2 

-1.4 

1.2 

10 

90 

'-0.7 

-0.1 

-1.4 

1.3 

15 

90 

-0.7 

0.0 

-1.3 

1.3 

20 

90 

-0.7 

0.1 

-1.4 

1.5 

25 

87 

-0.7 

0.0 

-1.6 

1.6 

30 

82 

-0.6 

0.0 

-1.3 

1.3 

35 

73 

-0.5 

0.0 

-1.2 

1.2 

40 

30 

-0.3 

0.6 

-1.1 

1.7 

45 

12 

-0.4 

0.0 

-1.0 

1.0 

50 

For  each  depth  of  a  bar  check  there  were  two  corrections  listed; 
the  first  correction  was  recorded  as  the  bar  was  being  lowered,  and  the 
second  during  raising.  These  two  corrections  were  averaged  to  obtain 
a  single  value  which  was  used  for  a  correlation  analysis  with  predicted 
tide  height. 
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A  plot  of  the  bar  check  and  tide  height  data  (Figure  2.4)  was 
made  and  visually  examined.  The  plot  indicates  a  poor  correlation 
between  bar  check  corrections  and  predicted  tide  heights.  For 
verification,  a  correlation  analysis  (Table  X)  was  performed  using 
methods  previously  discussed. 


TABLE  X 

Bar  Check  Statistical  Data  for  Survey  H-8422 


MEAN  TIDE  HEIGHT  .  3 . 3  f t 

STANDARD  DEVIATION  OF  TIDE  HEIGHT  .  2 . 9  f t  ' 

MEAN  CORRECTION  .  -0.7  ft  * 

STANDARD  DEVIATION  OF  CORRECTION  .  0.3  ft 

CORRELATION  COEFFICIENT  . ’ .  0.12  ! 

DEPTH  OF  OBSERVATIONS  .  25.0  ft 

t 


Because  of  the  poor  correlation,  actual  tidal  data  were  obtained  to 
investigate  the  possiblilty  that  the  results  were  related  to  the  predicted 
tides.  No  improvement  in  the  results  was  accomplished  by  using  actual 
tide  data.  No  further  manipulation  of  the  NOS  bar  check  data  was 
carried  out  because  of  the  number  of  unknowns  associated  with  the 
acquisition  procedures  and  because  the  amount  of  speculation  necessary 
could  not  lead  to  meaningful  analysis. 
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III.  DISCUSSION  OF  RESULTS 


A.  CREDDP  DATA 

Data  acquisition  procedures  of  CREDDP  made  an  investigation  of 
mean  sound  velocity  variation  possible,  but  there  are  factors  which 
make  the  study  incomplete.  Because  of  this,  a  complete  model  for  mean 
sound  velocity  prediction  and  correction  based  on  methods  presented  in 
this  thesis  was  not  fully  developed.  However,  the  primary  components 
of  such  a  model  were  devised  and  will  be  explained  in  this  chapter. 

The  factors  which  make  the  investigation  incomplete  and 
reopmmendations  concerning  verification  of  ideas  and  methods  presented 
in  this  thesis  will  be  addressed  in  Chapter  IV. 

Investigation  of  the  CREDDP  data  revealed  that  sound  velocity  in 
the  Columbia  River  estuai  y  varies  with  position,  depth,  and  time.  To 
analyze  the  problem  the  study  area  was  divided  into  two  primary 
sections,  the  north  channel  and  the  south  channel  (Table  XI).  The 
'river  mile'  distance  represents  the  shortest  deep-water  path  from  buoy 
"7"  at  the  river  entrance  to  the  station.  The  south  channel,  which  is 
maintained  for  commercial  marine  traffic,  has  a  river  mile  scale  which 
can  be  found  on  NOS  chart  18521,  whereas  north  channel  stations  were 
scaled  using  chart  18521. 

Two  categories  of  mean  sound  velocity  variation  are  considered. 
The  term  "horizontal  variation"  designates  differences  in  mean  sound 
velocity  between  stations;  i.e.,  how  the  velocity  changes  with  location. 
The  horizontal  variation  analysis  is  needed  to  develop  zones  or  regions 
where  velocity  correctors  may  be  grouped.  The  term  "temporal 
variation"  designates  the  change  in  mean  sound  velocity  with  time  at 
each  station.  With  temporal  variation,  mean  sound  velocity  correctors 
for  each  zone  change  with  time. 
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TABLE 

XI 

Channel 

Data 

north 

channel 

south 

channel 

Station 

River  Mile 

Station 

River  Mile 

2N 

2.0 

2S 

2.0 

4NB 

7.5 

4SB 

8.0 

5NB 

9.5 

5SA 

15.5  1 

1 

1 

j 

1 .  Horizontal  Variation 

For  the  analysis  of  horizontal  variatic^n  the  two  requirements 
made  were:  (1)  that  for  each  depth  there  be  but  one  horizontal 
variation  analysis  and;  (2)  the  average  of  all  mean  velocities  obtained  at 
the  depth  being  analyzed  be  designated  as  the  average  velocity.  Use 
of  the  average  velocity  allows  each  station  to  have  a  single  velocity  for 
each  depth.  The  average  velocity  presents  a  problem  if  more 
observations  were  made  on  one  end  of  the  tidal  spectrum  than  the 
other.  For  example,  if  more  casts  were  acquired  during  periods  of 
high  water  than  periods  of  low  water,  the  average  velocity  would  be 
higher  than  if  the  majority  of  casts  were  acquired  during  low  water 
periods.  One  of  the  CREDDP  objectives  was  to  acquire  casts  on  an 
hourly  basis.  This  objective  was  met  in  most  cases,  but  there  were 
some  instances  where  the  time  interval  between  successive  casts  was 
greater  than  1  hour.  For  example,  station  4SB  was  chosen  ^  'er  station 
4SA  for  designation  of  the  south  channel  because  there  were  no  casts 
acquired  at  station  4SA  during  the  higher  high  tide.  Another  potential 
problem  would  be  introduced  if  there  were  significant  variations  in  the 
river  flow.  According  to  Jay  [Ref.  5:  p.  114],  river  flow  remained  low 
for  the  acquisition  period.  An  example  of  horizontal  variation  of  the 
average  velocity  is  shown  in  Figure  3.1.  These  plots  indicate  the  rate 
of  horizontal  variation  of  the  average  velocity  is  approximately  the  same 


for  the  north  and  south  channels.  Data  for  the  11-  and  13-meter 
depths  were  limited  and  did  not  cover  the  range  of  tide  well  in  many 
cases. 

Regional  zoning  of  mean  sound  velocity  correctors  based  solely 
on  horizontal  variation  can  be  accomplished  with  little  difficulty.  Using 
the  depth  with  maximum  horizontal  variation,  zones  can  be  drawn  based 
on  the  NOS  criteria  of  plus  or  minus  4  meters  per  second.  A  major 
problem  still  exists  in  that  the  temporal  variation  of  mean  sound  velocity 
has  not  been  accounted  for  when  only  the  horizontal  variation  of  the 
average  velocity  is  considered. 

2.  Temporal  Variation 

An  analysis  of  temporal  variation  of  the  mean  sound  velocity 
demonstrated  that  proper  correction ^of  soundings  in  the  Columbia  River 
is  a  multi-dimensional  problem.  Sound  velocity  variables  change 
significantly  with  time  at  the  same  location,  and  the  magnitude  of 
variation  can  depend  on  a  number  of  cyclic  and  noncyclic  phenomena. 
Sound  velocity  variables  change  with  the  tide,  and  the  magnitude  of  the 
variation  is  affected  by  monthly  tidal  events  such  as  spring  and  neap 
tides.  The  temporal  and  horizontal  variations  of  velocity  for  the 
stations  of  the  north  and  south  channels  are  shown  in  Figure  3.1. 

Data  for  stations  2N  and  2S  were  acquired  during  a  neap  tide 
period  when  the  predicted  tide  range  was  minimum,  whereas  data 
acquisition  for  stations  4SA,  4SB,  5SA,and  5SB  was  accomplished 
during  a  spring  tide  period.  Data  for  stations  4NA,  4NB,  5NA,  5NB, 
and  5NC  were  acquired  during  a  period  of  average  tide  range.  As 
illustrated  by  Figure  3.1,  stations  4NB  and  4SB  have  about  the  same 
average  velocity  at  the  standard  depths,  but  with  respect  to  temporal 
variation,  the  range  of  velocity  for  station  4SB  was  significantly  larger 
than  that  found  at  station  4NB.  According  to  Jay  [Ref.  5;  pp 
114,117],  the  neap  tide  period  was  a  one  of  minimum  vertical  mixing  and 
maximum  stratification,  whereas  during  the  spring  tide  period,  maximum 
mixing,  minimum  stratification,  and  maximum  salinity  intrusion  occurred. 
Examination  of  the  casts  for  stations  4NB  and  4SB  (Table  XII) 
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SOUND  VELOCITY  RANGE 


O  AVG  VELOary-NORTVI  CHANNEL 
A  AVG  VEUJCITY-SOUTH  CHANh£L 

X  vELoarr  range 


4NB  5NB 

!S  4SB  5SA 

- 1 - 1 - 1 - T - 1 - 1 - 1 

0  4.0  6.0  8.0  10.0  12.0  14.0  16.0 

RIVER  MILE 


Temporal  and  Horizontal  Variations  of  Velocity 
Depth  =  9  m. 


confirmed  the  previous  statement.  The  cast  for  station  4NB  acquired  at 
predicted  low  tide  showed  a  significant  difference  in  temperature  when 
compared  to  the  low  tide  cast  of  station  4SB.  The  warmer  temperatures 
of  the  low  tide  cast  caused  the  the  minimum  velocity  to  be  higher  at 
station  4NB  than  at  station  4SB.  Comparison  of  the  high  tide  casts 
indicated  the  surface  layer  salinity  was  6.82  parts  per  thousand  less  at 
station  4NB  than  at  station  4SB.  The  relatively  large  salinity 
difference  for  the  high  tide  casts  caused  the  maximum  mean  velocity  to 
be  lower  at  station  4NB  when  compared  to  the  maximum  mean  velocity  at 
station  4SB.  Therefore,  the  increased  energy  associated  with  a  larger 
tide  range  can  increase  the  mean  sound  velocity  range  for  certain  areas 
in  two  ways.  First,  increased  mixing  associated  with  the  larger  tide 
range  can  moderate  temperature  for  low  tide  cast,  thereby  causing 
minimum  mean  velocities  to  be  lower.  Second,  increased  salinity 
intrusion  associated  with  a  larger  tide  range  can  cause  maximum  mean 
velocities  to  be  higher  when  compared  to  maximum  mean  velocities 
observed  during  an  average  or  neap  tide  range. 

a.  Correlation  of  Tide  Height  and  Sound  Velocity 

The  effort  to  correct  mean  sound  velocity  for  temporal 
variation  began  with  a  correlation  analysis  of  predicted  tide  height  and 
mean  sound  velocity.  The  results  (Table  XIII)  indicate  a  strong 
relationship  between  tide  height  and  mean  sound  velocity  for  most  of 
the  CREDDP  data.  Correlation  in  the  1 -meter  layer  was  usually  not  as 
strong  as  the  deeper  layers  for  most  of  the  stations.  Dashed  line 
entries  in  Table  XIII  indicate  insufficient  data  for  a  correlation 
analysis . 

b.  Regression  Analysis 

First-  (linear),  second-,  and  third-order  polynomial 
regression  analyses  were  performed  for  each  standard  depth  of  the 
CREDDP  sampling  stations.  The  objective  of  these  analyses  was  to 
provide  a  numerical  method  of  predicting  mean  sound  velocity  using  the 
predicted  tide  height.  These  methods  would  be  considered  completely 


TABLE  XII 

The  Effect  of  Tide  Range  on  Mean  Sound  Velocity  Range 


STATION  4NB 


TIDE  HEIGHT  =  0. 

02  ft 

(LOW  TIDE) 

DEPTH 

TEMP 

SAL 

VEL 

m 

g/kg 

m/s 

1.0 

14.22 

6.32 

1472.2 

3.0 

14.  14 

7.48 

1472.5 

5.0 

14.05 

8.52 

1473 . 1 

7.0 

13.93 

9.73 

1473.5 

TIDE  HEIGHT  =  8. 

3  ft  (HIGH  TIDE) 

DEPTH 

TEMP 

SAL 

VEL 

m 

g/kg 

m/s 

1.0 

12.18 

23.49 

1484.6 

3.0 

11.78 

29.26 

1487.4 

5.0 

11.97 

29.89 

1489.2 

7.0 

11.66 

30.19 

1489.8 

STATION  4SB 

TIDE  HEIGHT  =  -1 

,.5  ft 

(LOW  TIDE) 

DEPTH 

TEMP 

SAL 

VEL 

m 

g/kg 

m/s 

1.0 

12.71 

8.41 

1468.6 

3.0 

12.69 

8.92 

1468.9 

5.0 

12.66 

9.55 

1469.2 

7.0 

12.64 

9.70 

1469.4 

TIDE  HEIGHT  =  9. 

2  ft  (HIGH  TIDE) 

DEPTH 

TEMP 

SAL 

VEL 

m 

g/kg 

m/s 

1.3 

11.88 

30.31 

1491.9 

3.0 

11.87 

30.69 

1492 . 1 

5.0 

11.84 

31.35 

1492.5 

7.0 

11.84 

31.48 

1492.7 

TABLE  XIII 

Correlation  Coefficients 


STN 

Im 

3m 

5m 

DEPTH 

7m 

9m 

11m 

13m 

2N 

0.61 

0.81 

0.85 

0.87 

0.87 

0.87 

0.78 

2S 

0.86 

0.68 

0.79 

0.83 

0.88 

0.89 

0.85 

4NA 

0.83 

0.91 

0.94 

0.97 

0.98 

0.99 

0.99 

4NB 

0.81 

0.90 

0.94 

0.96 

0.96 

0.99 

.... 

4SA 

0.71 

0.82 

0.86 

0.89 

0.92 

0.90 

0.93 

4SB 

0.75 

0.83 

0.89 

0.93 

0.96 

0.99 

.... 

SNA 

0.86 

0.89 

0.92 

0.91 

0.90 

_  ■  -  - 

.... 

5NB 

0.80 

0.86 

0.91 

0.94 

0.95 

0.97 

0.96 

5NC 

0.90 

0.91 

0.91 

0.94 

0.96 

0.97 

0.97 

5SA 

0.89 

0.88 

0.88 

0.87 

0.85 

0.86 

.... 

5SB 

0.88 

0.88 

0.88 

0.87 

0.86 

-  -  _  - 

.... 

regression  analysis  was  for  comparison  of  fit,  as  well  as  comparison  of 
the  predicted  velocities  and  associated  confidence  intervals. 

A  potential  source  of  error  may  exist  when  using 
predictions  calculated  from  tide  height  values  which  are  outside  of  the 
range  of  values  used  for  the  initial  regression  analysis.  For  example, 
if  all  casts  were  acquired  during  a  neap  tide  period,  such  as  was  the 
case  for  stations  2N  and  2S,  velocities  predicted  by  the  regression 
analysis  for  high  or  low  tide  heights  of  a  spring  tide  period  are  subject 
to  error.  For  stations  4SB  and  4SA,  predictions  were  calculated  using 
a  9.0-foot  tide  height;  station  4SB  (Figure  3.2)  was  occupied  during 
higher  high  tide,  whereas  station  4SA  (Figure  3.3)  was  not  occupied 
during  this  period.  The  second-  and  third-order  analysis  predicted 
velocities  which  probably  were  not  physically  possible  at  station  4SA. 
Therefore,  a  basic  objective  in  using  polynomial  regression  models  for 
mean  sound  velocity  predictions  and  corrections  should  be  the 
determination  of  the  maximum  and  minimum  velocities  expected,  and  the 
minimum  range  of  tide  where  these  velocities  could  occur.  The  most 
accurate  method  for  identification  of  these  values  would  be  to  measure 
the  sound  velocity  variables  at  the  same  location  covering  more  phases 


of  the  lunar  cycle.  The  danger  of  predicting  outside  of  the  observed 
range  is  illustrated  by  Figure  3.3,  whereas  Figure  3.2  shows  that 
predicting  within  the  observed  range  can  be  done  accurately. 

Determination  of  the  most  efficient  polynomial  regression 
analysis  to  use  for  prediction  of  mean  sound  velocity  was  accomplished 
by  examination  of  plots  of  the  data  arrays  and  comparison  of  the 
calculation  of  predicted  velocities  and  confidence  intervals  provided  by 
each  analysis. 

(1)  Data  Array  Plots.  Examination  of  the  data  array 
plots  indicated  that  the  first-order  (linear)  analysis  would  be  sufficient 
for  most  of  the  data  arrays,  although  some  of  the  arrays  did  suggest  a 
curvilinear  trend.  Curvature  could  be  detected  in  the  plots  of  stations 
4NA,  4SB  (Figure  3.2),  4SA  (Figure  3.3),  and  4NB  (Figure  3.4),  with 
maximum  curvature  in  the  plots  of  stations  4SA  and  4SB.  This  may  be 
an  indication  that  the  degree  of  curvature  is  proportional  to  the  range 
of  tide  since  data  for  stations  4SA  and  4SB  were  acquired  during  the 
spring  tide  period,  whereas  data  acquisition  for  stations  4NA  and  4NB 
was  carried  out  during  a  period  of  average  tide  range.  This 
hypothesis  cannot  be  verified  without  further  investigation. 

(2)  Predicted  Velocities .  A  computer  program  was 
written  to  predict  the  mean  velocity  for  five  different  heights  of  tide 
with  the  maximum  and  minimum  heights  being  the  same  as  those  used  in 
the  regression  analysis.  The  other  three  heights  were  equally  spaced 
between  the  maximum  and  minimum.  Rarely  was  there  a  significant 
difference  (plus  or  minus  4  meters  per  second)  between  the  predictions 
calculated  by  the  three  orders  of  polynomial  regression  (Table  XIV). 
Mean  velocities  predicted  using  the  the  third-order  regression  analysis 
were  in  error  for  stations  2S  and  5NB  at  the  13-meter  depth.  The 
cause  of  this  error  could  not  be  explained. 

(3)  Confidence  Intervals .  Confidence  intervals  (Table 
XV)  were  computed  for  the  predicted  velocities  using  t  values  for 
90-percent  confidence.  In  general  terms,  the  confidence  intervals  were 
larger  for  all  stations  in  the  1-  and  3-meter  depths,  with  the  intervals 
becoming  less  with  increasing  depth.  This  corresponds  with  the  lower 


TABLE  XIV 

Predicted  Velocities- -Depth  =  5  m 


STN 

2N 

1.2 

TIDE 

2.6 

HEIGHT  (ft) 

4.0  5.5 

6.9 

ORDER 

1 

2 

3 

1475.5 

1474.8 

1474.1 

1479.0 

1478.9 

1479.3 

1482.5  . 

1482.7 

1482.5 

1486.2 

1486.2 

1485.3 

1489.7 
1489.2 

1488.7 

STN 

2S 

0.3 

TIDE 

1.9 

HEIGHT  (ft) 

3.6  5.2 

6.9 

ORDER 

1 

2 

3 

1476.8 

1479.5 

1479.3 

A 

1480.0 

1479.6 

1479.8 

1483.3 

1481.9 

1481.9 

1486.5 

1486.1 

1486.0 

1489.9 
1492.8 
1493 . 1 

STN 

4NB 

0.0 

TIDE 

2.1 

HEIGHT  (ft) 

4.1  6.2 

8.3 

ORDER 

1 

2 

3 

1472.8 

1474.4 

1474.4 

1476.6 

1475.9 

1476.0 

1480.3 

1478.8 

1478.8 

1484.2 

1483.5 

1482.9 

1488.0 

1489.6 

1488.5 

STN 

4SB 

-1.5 

TIDE 

1.2 

HEIGHT  (ft) 

3.8  6.5 

9.2 

ORDER 

1 

2 

3 

1465.0 

1469.2 

1469.5 

1471.5 

1471.1 

1470.9 

1477.8 

1475.9 
1475.6 

1484.3 

1483.9 

1482.9 

1490.8 

1494.9 
1491.4 

STN 

5NB 

0.2 

TIDE 

2.2 

HEIGHT  (ft) 

4.2  6.2 

8.3 

ORDER 

1 

2 

3 

1468.9 

1470.5 

1471.1 

1472.4 
1472.0 

1471.5 

1476.0 

1475.0 

1475.3 

1479.6 
1479.3 

1480.6 

1483 . 4 
1485.2 
1485.7 

STN 

5SA 

-1.6 

TIDE 

1.2 

HEIGHT  (ft) 

4.0  6.7 

9.5 

ORDER 

1 

2 

3 

1460.2 
1461 . 6 
1461.9 

1463.1 

1462.7 

1462.5 

1465.9 

1465.0 

1464.7 

1468.6 

1468.3 

1467.2 

1471.5 
1473 . 1 
1468.8 

ORDER  =  Degree  of  polynomial  regression 


correlation  coefficients  found  in  the  shallow  depths.  Larger  confidence 
intervals  were  also  found  for  data  acquired  at  stations  2N,  2S,  4SA, 
and  4SB.  These  stations  had  lower  correlation  coefficients  and  higher 
residuals  when  compared  to  the  other  stations.  For  stations  2N  and  2S, 
this  may  indicate  less  stability  due  to  the  close  proximity  of  the  open 
ocean.  As  stated  before,  data  for  stations  4SA  and  4SB  were  acquired 
during  spring  tide  conditions,  which  tends  to  indicate  less  stability  for 
prediction. 

In  most  cases,  increasing  the  degree  of  the 
regression  analysis  did  not  produce  more  accurate  results  in  terms  of 
confidence  intervals.  The  standard  error  of  estimate  was  sometimes 
lower  with  a  higher  order  analysis,  but  the  reduction  in  the  degrees  of 
freedom  countered  the  effect  of  the  lower  standard  error  of  estimate  for 
computatipn  of  confidence  intervals. 

Confidence  interval  results  indicate  that  correcting 
soundings  for  mean  velocity  variation  based  on  predictions  provided  by 
the  regression  analysis  cannot  be  considered  a  complete  success. 
Attaining  the  goal  of  predictions  to  within  plus  or  minus  4  meters  per 
second  may  not  be  possible  in  some  cases,  but  the  regression  analysis 
does  show  a  significant  reduction  in  error  when  compared  to  only  using 
average  velocity  and  horizontal  variation  as  the  correction  procedure. 

B.  NOS  BAR  CHECK  DATA 

Analysis  of  the  bar  check  data  revealed  too  many  inconsistencies 
and  required  too  much  speculation  for- comparison  of  the  CREDDP  data 
results.  A  major  problem  of  the  direct  comparison  method  of  the  bar 
check  is  that  examination  of  the  sound  velocity  variables  and  any 
change  in  their  magnitude  is  not  possible.  The  bar  check  data  did 
indicate  the  possibility  of  either  an  extremely  high  sound  velocity 
variation  or  an  instrument  error  which  was  not  constant  with  depth. 

Most  likely,  the  results  show  an  undeterminable  combination  of  sound 
velocity  variation  and  instrument  error.  The  problem  of  using  the  bar 
check  data  for  any  meaningful  analysis  is  illustrated  by  Table  XVI. 


TABLE  XV 


Confidence  Intervals --Depth  =  5  m 


STN 

2N 

1.2 

TIDE 

2.6 

HEIGHT  (FT) 
4.0 

5.5 

6.9 

ORDER 

1 

4.2 

3.9 

3.8 

3.8 

4.0 

2 

4.  1 

3.9 

3.7 

3.7 

3.9 

3 

4.3 

4.0 

3.9 

3.9 

4.1 

STN 

2S 

0.3 

TIDE 

1.9 

HEIGHT  (FT) 
3.6 

5.2 

6.9 

ORDER 

1 

4.3 

4.1 

4.0 

4.1 

4.3 

2 

3.8 

3.6 

3.5 

3 . 6 

3.7 

3 

3.9 

3.7 

3 . 6 

3.7 

3.8 

STN 

4NB 

0.0 

TIDE 

2.1 

HEIGHT  (FT) 
4.1 

6.2 

8.3 

ORDER 

1 

3.1 

2.9 

2.8 

2.8 

3.0 

2 

2.9 

2.7 

2.6 

■  2.7 

2.8 

3 

3.3 

3.1 

3.0 

3.0 

3.2 

STN 

4SB 

-1.5 

TIDE 

1.2 

HEIGHT  (FT) 

3 . 8 

6.5 

9.2 

ORDER 

1 

5.2 

5.0 

4.9 

4.9 

5.1 

2 

4.9 

4.7 

4.6 

4.6 

4.8 

3 

4.9 

4.7 

4.6 

4.6 

4.8 

STN 

5NB 

0.2 

TIDE 

2.2 

HEIGHT  (FT) 
4.2 

6.2 

8.3 

ORDER 

1 

3 . 3 

3.2 

3.2 

3.2 

3.3 

2 

3.0 

3.0 

2 . 9 

2 . 9 

3 . 1 

3 

2.9 

2.8 

2.8 

2.8 

2.9 

STN 

5SA 

-1.6 

TIDE 

1.2 

HEIGHT  (FT) 
4.0 

6.7 

9.5 

ORDER 

1 

3.9 

3.8 

3.7 

3.8 

3.9 

2 

3.6 

3.5 

3.5 

3.5 

3 . 6 

3 

4.5 

4.4 

4.4 

4.4 

4.6 

units  =  ±  tn/s 


TABLE  XVI 

Bar  Check  Correction  Data  for  Survey  H-8422 


UPPER  LOWER 


DBS 

MEAN 

ft 

LIMIT 

ft 

LIMIT 

ft 

RANGE 

ft 

DEPTH 

ft 

DEL  V 
m/s 

18 

-1.0 

-0.6 

-1.4 

0.8 

5 

231.5 

0 

-0.8 

-0.2 

-1.4 

1.2 

10 

177.5 

'0 

-0.7 

-0.1 

-1.4 

1.3 

15 

124.4 

10 

-0.7 

0.0 

-1.3 

1.3 

20 

93.5 

0 

-0.7 

0.1 

-1.4 

1.5 

25 

88.3 

17 

-0.7 

0.0 

-1.6 

1.6 

30 

78.4 

12 

-0.6 

0.0 

-1.3 

1.3 

35 

53.5 

'3 

-0.5 

0.0 

-1.2 

1.2 

40 

44 . 4 

0 

-0.3 

0.6 

-1.1 

1.7 

45 

55.5 

DEL  V  represents  the  magnitude  of  velocity  difference 
necessary  to  compensate  for  the  range  of  the 
corrections  if  there  were  no  instrument  error. 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 


Investigation  of  sound  velocity  variation  in  the  Columbia  River 
estuary  revealed  the  following  significant  findings: 


Sound  velocity  variation  in  certain  areas  required  the  examination 
of  how  veloctiy  changed  with  time  and  location  if  NOS  accuracy 
standards  were  to  be  met.  The  procedure  of  computing  only  one 
corrector  table  based  on  an  average  velocity,  or  correction,  was 
not  acceptable. 


CREDDP  data  analysis  showed  a  strong  relationship  between  tide 
height  and  sound  velocity,  which  stimulated  the  development  of  a 
method  to  account  for  the  variation  of  sound  velocity  with  time. 
Using  a  regression  analysis,  the  accuracy  of  velocity  correctors 
was  improved  significantly,  although  there  were  still  instances 
where  NOS  standards  could  not  be  met. 


NOS  bar  check  data  analysis  showed  a  significant  variation  in 
corrections,  although  the  variation  did  not  correlate  with  tide 
height  as  the  CREDDP  data  did.  Analysis  of  the  variation  found 
at  each  5-foot  increment  of  depth  revealed  inconsistencies  which 
could  not  be  explained. 


B.  RECOMMENDATIONS 


The  following  recommendations  are  made: 


Further  study  of  sound  velocity  variation  must  be  carried  out  for 
verification  of  ideas  and  methods  presented  in  this  thesis. 

Analysis  of  data  acquired  over  longer  time  periods  at  the  same 
locations  will  allow  a  more  complete  study  of  the  variation 
associated  with  different  lunar  phases  and  the  effect  of  changes  in 
river  flow.  NOS  CTD  data  were  acquired  in  the  Columbia  River 
estuary  during  the  same  time  period  as  CREDDP;  these  data  may 
allow  study  or  the  sound  velocity  variation  over  longer  time 
periods . 


NOS  bar  check  data  analysis  showed  that  the  hydrographer  should 
be  equipped  with  an  efficient  method  of  measuring  sound  velocity 
in  addition  to  determining  instrument  error.  Knowing  the  value  of 
sound  velocity,  detection  of  significant  changes  and  trends  is  more 
efficient  and  accurate.  Bar  checks  should  not  be  used  as  the 
primary  method  of  correcting  soundings,  but  should  be  reserved 
for  determination  of  instrument  error. 


Hydrographers  conducting  surveys  in  areas  where  reversing  tides 
could  cause  a,  significant  variation  in  sound  velocity  may  benefit  by 
noting  the  height  of  tide  for  each  cast  or  bar  check  acquired  for 
determining  sound  yelocitv  corrections.  If  significant  variation  is 
detected,  the  possibilty  of  predicting  and  correcting  for  temporal 
variation  .using  the  tide  height  as  an  independent  variable  should 
be  investigated.  FORTRAN  computer  programs  written  for  the 
methods  presented  in  this  thesis  are  listed  in  Appendices  A,  B,  C, 
D,  E,  and  F. 


In  areas  where  temporal  variation  of  sound  velocity  is  significant, 
the  hydrographer  may  benefit  by  using  a  sound  velocimeter  to 
correct  soundings.  The  velocimeter  would  require  only  a  few 
minutes  to  deploy  and  retrieve.  The  hydrographer  could  determine 
sound  velocity  periodically  throughout  the  working  day  without 
significantly  hindering  operations. 


APPENDIX  A 

FORTRAN  PROGRAM  PACKAGE- -  SOUND  VELOCITY  COMPUTATIONS 
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DO  50  K=2,J 

SUM=SUM+((D(K)-D(K-1))*(SVL(K)+SVL(K-1)))/(2*SVL(K)*SVL(K-1)) 

CONTINUE 

HVL(J)=(D(J)-D(1))/SUM 
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FORTRAN  PROGRAM  PACKAGE- -BASIC  STATISTICS 
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VELOCITY  IS  ALSO  COMPUTED.  A  HIGH  CORRELATION  COEFFICIENT  (R  > 
MAY  INDICATE  THAT  SOUND  VELOCITY  CAN  BE  PREDICTED  BASED  ON  TIDE 
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R=  (  (N-'SUMVT)-  (SUMT*SUMV)  )/SQRT(  (  (N*SUMT2  )  -  ( SUMT**2  )  )*(  (N*SUMV2)- 
*(SUMV**2))) 

PRINT  STATION,  AVERAGES,  STANDARD  DEVIATIONS,  AND  CORRELATION 
COEFFICIENT. 

WRITE (7,110)  STN ( J ) , AVT , SDTI , AVV , SDVI , R 


110  FORMAT ( IX , A3 , 2X . F5 . 1 , 2X , F5 . 1 , 2X , F7 . 1 , 2X , F5 . 1 , 2X , F5 . 2 ) 

5  CONTINUE 


FORTRAN  PROGRAM  PACKAGE- -FIRST-ORDER  REGRESSION. 
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SOUND  VELOCITY  AND  X  IS  THE  TIDE  HEIGHT.  IF  THERE  IS  GOOD 
CORRELATION  BETWEEN  TIDE  HEIGHT  AND  SOUND  VELOCITY,  THIS  PROGRAM 
MAY  ALLOW  HYDROGRAPHERS  TO  CORRECT  FOR  THE  TEMPORAL  VARIATION  OF 
SOUND  VELOCITY . 
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SUM2  =  SUM2  + ( ER ( I ) **2 ) 

PRINT  OBSERVATION  NUMBER  AND  RESIDUAL. 
WRITE(7,110)  I,ER(I) 

110  F0RMAT(16X,I2,11X,F6.2) 

30  CONTINUE 


PRINT  SUM  OF  RESIDUALS,  SUM  OF  RESIDUALS  SQUARED,  AND  NUMBER  OF  OBS. 
WRITE (7, 125)  SUM1,SUM2,B 
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A3=0.0 

READ  IN  TIDE  HEIGHTS  AND  OBSERVED  SOUND  VELOCITIES.  A  FLAG  OF 
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A2=(G-((F*J)/H)^(A1*(((D*F)/H)-C)))/L 
A3= ( J- (A1*D) - (A2*F) )/H 

PRINT  STATION  ID  AND  COEFFICIENTS  OF  PREDICTION  EQUATION. 
V/RITE(7,120)  STN,A1,A2,A3 


120  FORMAT(7X,A3,3X,3F7.2) 
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APPENDIX  F 

FORTRAN  PROGRAM  PACKAGE- -PREDICTED  VELOCITY/ CONFIDENCE  INTERVAL 
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